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Abstract 

We extensively develop a perturbation theory for nonlinear cosmological dynamics, based on the 
Lagrangian description of hydrodynamics. We solve hydrodynamic equations for a self-gravitating 
fluid with pressure, given by a polytropic equation of state, using a perturbation method up to 
second order. This perturbative approach is an extension of the usual Lagrangian perturbation 
theory for a pressureless fluid, in view of inclusion of the pressure effect, which should be taken 
into account on the occurrence of velocity dispersion. We obtain the first-order solutions in generic 
background universes and the second-order solutions in wider range of a polytropic index, whereas 
our previous work gives the first-order solutions only in the Einstein-de Sitter background and the 
second-order solutions for the polytropic index 4/3. Using the perturbation solutions, we present 
illustrative examples of our formulation in one- and two-dimensional systems, and discuss how the 
evolution of inhomogeneities changes for the variation of the polytropic index. 
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I. INTRODUCTION 



Hydrodynamics is a powerful tool to study various astrophysical phenomena, from those 
associated with compact objects up to large-scale structure formation. For example, when 
investigating gravitational instability of cold dark matter for structure formation, analyses 
are easier to handle by adopting a hydrodynamical description, rather than by trying to solve 
the Boltzmann equation of a self-gravitating N-particles system. The linear perturbation 
theory of a homogeneous and isotropic universe ]T|, 0, [|, [J is a typical case, which gives 
a qualitative estimate for gravitational instability. It is based on the Eulerian picture of 
hydrodynamics, while approximations based on the Lagrangian hydrodynamics has been 
recognized to be more useful, such as the celebrated Zel'dovich approximation || |7|, ||. 
This paper deals with an approximation theory of gravitational instability based on the 
Lagrangian hydrodynamics. 

Although the Zel'dovich approximation has been found to give an accurate description up 
to the stage where density perturbations grow to be unity, it involves a serious shortcoming 
that it cannot be applied after caustics in the density field are formed. In the Zel'dovich 
approximation, the fluid elements continue to move in the directions that are determined 
by initial conditions all the time, and consequently high density regions such as 'pancakes' 
cannot stay compact beyond the caustic formation, while numerical simulations have shown 
the presence of clumps with a very wide range in mass at any given time || . Moreover, once 
caustics in the density field are formed, a hydrodynamical description itself is not valid in 
general. Then, do we have to abandon a hydrodynamical description and try to solve the 
Boltzmann equation, or tackle N-body simulations? 

In order to proceed with a hydrodynamical description without the formation of caustics, 
the 'adhesion approximation' |l(J has been proposed, in which an artificial viscosity term is 
added to the Zel'dovich approximation. This modified approximation successfully describes 
the stage where the original Zel'dovich approximation breaks down, but the physical origin 
of the viscosity term should be clarified. Some remarkable works have been done on this 
issue; Buchert and Dommguez [O] argued, by beginning with the collisionless Boltzmann 
equation |12|], that the effect of velocity dispersion become important beyond the caustics, 
where multi-stream flow arises, and the presence of velocity dispersion of the flow can yield 



pressure-like or viscosity terms; Buchert et al. [13] showed how the viscosity term is generated 
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by a pressure-like force of a fluid under the assumptions that the peculiar acceleration is 



parallel to the peculiar velocity; Dommguez [TJj clarified that a hydrodynamic formulation is 
achieved via a spatial coarse graining in a many-body gravitating system, and the 'adhesion 
approximation' can be derived by the expansion of coarse-grained equations with respect to 
the smoothing length. (See also Ref. jn|.) In these works, they also obtained implications for 
an 'equation of state', which is a phenomenological relationship between kinematical pressure 
P and mass density p. Buchert and Dommguez [|TTJ found that, if the effect of velocity 
dispersion is small and the velocity dispersion is approximately isotropic, the equation of 
state should take the form P oc p 5//3 ; Buchert et al. [|l3j showed that an adhesion- like equation 
can be derived if the equation of state is assumed as P oc p 2 ; Moreover, a plausible value of 
the polytropic index 7 = dlnP/dlnp has been found to be close to 5/3 from cosmo logical 



N-body simulations [16]. 

From these aspects, it is of interest to extend a Lagrangian perturbation scheme to a fluid 
with pressure, and to explore how the scheme works as an approximation for cosmological 
structure formation. Actually, Adler and Buchert [[17) and Morita and Tatekawa |L8| have 
formulated perturbation theory in the Lagrangian hydrodynamics, taking into account the 
pressure effect under the assumption that the pressure is a function of the mass density only. 
In our earlier work , imposing a polytropic relation as the equation of state, we solved 
the Lagrangian perturbation equations up to second order for cases where the equations 
are solved easily, and showed illustrative examples with the solutions in a one-dimensional 
system. In particular, the second-order solutions were obtained only for the case in which the 
polytropic index 7 is 4/3, while a plausible value of 7 seems to be larger, as was mentioned 
above. 

In this paper, we extend our earlier work by solving the first-order perturbation equations 
in generic background universes and the second-order perturbation equations for wider range 
of a polytropic index, and by presenting illustrations in one- and two-dimensional systems. 
We examine how the behavior of the perturbation solutions, and the resultant evolution 
of inhomogeneities, change for the variation of the polytropic index. This enables us to 
discuss whether, or for what kind of the equation of state, the adhesion-type approximation 
is realized in the Lagrangian perturbation scheme. 

This paper is organized as follows. In Sec. II, we present Lagrangian hydrodynamic equa- 
tions, governing the system we consider. In Sec. fTJ, the first-order perturbation equations 
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are derived and their solutions are shown, not only in the Einstein-de Sitter background but 
also in more generic backgrounds. In Sec. |TV| , we obtain the second-order perturbation equa- 
tions and present their solutions in an approximate form for 7 > 4/3. Section |V] provides 
illustrative examples of our formulation in one- and two-dimensional models. In Sec. |VT|, we 
discuss our results and state our conclusions. 



II. BASIC EQUATIONS 

In this section, we present hydrodynamic equations in the Lagrangian description, which 
our approach stands on. The matter model we consider is a self-gravitating fluid with energy 
density p and pressure P, which arises in the presence of velocity dispersion. Then the basic 
equations we start from are 

dp n a 1„ , s 

7 £ + 3-p + -V iE - (pv =0, 1 
at a a 

dv a 1 . „ . 1 „ „ 

1 - + -v + -(vV x )v = g V X P, (2) 

at a a pa 

V x xg = 0, (3) 

^ x g = -4nGa(p - p b ) , (4) 

where v and g are the peculiar velocity and the peculiar gravitational field, respectively, 
which represent the deviation from a background, homogeneous and isotropic universe. The 
cosmic scale factor a(t) and the energy density phit) of the background universe satisfy the 
Friedmann equations 

a\ 2 8nG K. A 

a) 3 a z 3 

a AttG A , s 

a = -— " b+ 3- (6) 

with a curvature constant K, and a cosmological constant A. In order to solve the hydrody- 
namic equations, we must specify an equation of state. Throughout this paper, we consider 
barotropic fluids, in which the pressure P is a function of the energy density only, P = P(p). 
Introducing the Lagrangian time derivative 



at at a 
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Eqs. (|1]) and (||) become 

^ + 3^+^.-10 = 0, (7) 
at a a 

dv a 1 „ „ 

— + -v = g V X P. 8 

at a pa 

In the Lagrangian hydrodynamics, the coordinates x of the fluid elements are represented 
in terms of Lagrangian coordinates q as 

x = q + s(q,t), (9) 

where q are defined as initial values of x, and s denotes the Lagrangian displacement vector 
due to the presence of inhomogeneities. The exact form of the energy density is then obtained 
from Eq. ([/]) as 

P = PbJ" 1 , (10) 

where J = det(dxi/dqj) = det(Sij+dsi/dqj) is the Jacobian of the coordinate transformation 
from x to q. The peculiar velocity is v = as, and from Eq. (U), the peculiar gravitational 
field is written as 

where an overdot (') denotes d/dt. Hence, from Eqs. and (H), we obtain the following 
equations for s: 

V x x (s + 2^sj =0, (12) 

V x -[s + 2-s-- — (p) J- l V x J = -AnGp^J- 1 - 1) . (13) 
V a a 1 dp / 



If we find solutions of Eqs. (|12|) and ( |T3| ) for s, the dynamics of the system considered is 
completely determined. Since these equations are highly nonlinear and hard to solve exactly, 
we will advance a perturbative approach. Remark that, in solving the equations for s in 
the Lagrangian coordinates q, the operator Va; will be transformed into V 9 by the following 
rule: 

d dxj d d dsj d . d . 

' ■ TT" 1 — • (14) 



dqi dqi dxj dxi dqi dx 



j 
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III. FIRST-ORDER SOLUTIONS 



Hereafter we develop a perturbative approach for the Lagrangian displacement vector s 
of the fluid elements. In the first-order approximation, Eqs. (12) and (|T3|) become 

Vq x (a^ + 2-a^)=0, (15) 



V, ■ (V) + - ~M V q (V q • = 47rG Pb V q • , (16) 

where denotes the first-order displacement vector in the perturbative expansion. De- 
composing into the longitudinal and the transverse modes as s^ 1 ' = V q S + S T with 



Vq • S T = 0, we have 



V q x (5 T + 2-5 T ) =0, (17) 



S + 2-S -AnGp h S - — — (p b ) V 2 q S = 0. (20) 



V 2 q (s + 2- a S - AnG Ph S - ^(Pb) V\s\ = . (18) 
These equations are reduced by imposing some adequate boundary conditions to 

S T + 2-S T = 0, (19) 
a 

l-S — AnGpbS ~ 

a a 2 dp 

In our previous paper, we obtained the perturbation solutions only for the Einstein-de 
Sitter background. Here we solve the equations for the first-order perturbations in generic 
background universes. Equation ( |T9"D can be integrated easily even in this case, although an 
explicit form of the solutions is not presented here. For Eq. (pO"D, the Fourier transformation 
with respect to the Lagrangian coordinates q yields 

§ + 2-§-4nG Ph S + \^(p h )\K\ 2 S = 0, (21) 
a a z dp 

where (■) denotes the Fourier transform, and K is a wavenumber associated with the La- 
grangian coordinates. Replacing the time variable t with a and using the Friedmann equa- 
tions (0) and @, we have 

87rGp b 2 „ A 2 \ d 2 S* / 2/C . \ dS ( 1 dP, , .„ I2 . „ \ & n 

— — ^a 2 - K + -a 2 ) —+ 4irG Ph a + Aa — + -j-(pb) \K\ 2 - AnGp h )S = 

6 6 J dcr V a / da \a z dp J 

(22) 
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If we assume a polytropic equation of state P = np" 1 with a constant k and a polytropic 
index 7, this equation becomes 

W ^A^ + m-^ + Aa )i + (^.^)s, , (23) 



V a 3 / da 2 V a 2 a / da y a 3 "' 1 a 3 

where Ci = A-kG p\,{a- in ) and C 2 = /«7Pb( a in) 7 ~ 1 'V*' 7 1 \ 

Let us consider solving Eq. (Effi). In the Einstein-de Sitter background, where /C = and 
A = 0, the solutions of Eq. (^) are written in a relatively simple manner. They are, for 
7 ^ 4/3, 

7 (4-3 7 )/2 



§(K,a) oc a 1/4 J± 5 /(8-6 7 ) fy^" 



2C 2 


K\ 






4 


-H 



(24) 



where J v denotes the Bessel function of order u, and for 7 = 4/3, 

S(K, a) oc fl -Vti\^WW . (25) 



In the nonflat backgrounds with /C 7^ and A = 0, the solutions of Eq. (|23| ) for 7 = 1, 4/3 
can be written in terms of Gauss' hypergeometric function T as 

S(K,a) oc a' 3 J 7 (a^, a 2 , a 3 ; 7^~r] > ( 26 ) 

where 



1 C 2 |K| 2 1 C 2 |K| 2 3 3 

(ai,a 2 ,a 3 ,/3) = ( -1 + W- + , -1 - d- + , -- , -- 



3 1 C 2 |K| 2 3 1 CJKl 2 7 . 



3 , 25 C 2 |K| 2 1 , 25 C 2 \K\ 



2 



(«!, a 2 , a 3 , /?) = - ± 4/ — — — , -- ± \ — 



4 V 16 2Ci 4 V 16 2Ci 
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25 2CW _1 25_0^1 
V 4 C x 4 V 16 2Ci 1 



In the flat (K = 0) backgrounds with A 7^ 0, we can also write the solutions of Eq. ( |2"3"D for 
7 = 1/3, 4/3 in the form 



S(K,a) (xaP T I ai,a 2 ,a 3 ; -^~r 



(29) 
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where 



(a 1 ,a 2 ,a 3 ,/3) = -J + VH 



1 /l CIKI 2 1 /l C 2 |K| 2 1 3 



6 V 9 3A'6 V 9 3A'6'2 



2 1 C 2 |K| 2 2 1 g 2 |JCp 11 

3 + V9-^A-'3"V9 SA-'Y' 1 ' fo ^ = 1 /3,(30) 



(ai,a 2 ,a 3 ,/3) 

9Ci 4 V 16 2d 

Let us note the relation between the behavior of the above solutions and the Jeans 
wavenumber, which is defined as 

1/2 




Ki = 



( AnGp h a 2 Y 
ldP/dp(p b )] 



The Jeans wavenumber, which gives a criterion whether a density perturbation with a 
wavenumber will grow or decay with oscillation, depends on time in general. If the polytropic 
equation of state P = up 1 is assumed, 



K, = ^ a^~ 4 )/ 2 . (32) 

Equation (JJ^) implies that, if 7 < 4/3, Kj will be infinitesimal and density perturbations 
with any wavenumber will decay in process of time, and if 7 > 4/3, all density perturbations 
will grow to collapse. This is confirmed by the form of the solutions, Eq. (0), by rewriting 
it as 

' ^ ^ (33) 



S{K,a) oc a 1/4 ^±5/(8-67) h 



4-3 7 | Kj J 

However, this fact seems to be curious because one may expect that, as the polytropic 
index 7 is larger, the effect of the pressure would be stronger and consequently the growth 
of density perturbations would be supressed more effectively. The unexpected result may 
be caused by construction of the first-order approximation, in which the strength of the 
pressure effect is determined only by the coefficient (1/a 2 ) dP/dp(pt>) in the fourth term of 
the left side of Eq. ( p0|) . The square of the 'sound speed,' dP/dp, which is contained in 
the coefficient, is originally a function of p, but now in the coefficient p is replaced with p b 
because of the first-order approximation. Since pb oc a -3 , the coefficient decays sooner as 
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the index 7 is larger, and it leads to the consequence. This problem may be resolved by 
trying higher-order approximations, where the pressure effect is provided not only by the 
background density but also by the presence of inhomogeneities. Let us proceed to second 
order, noticing the above fact. 

We should also note that the above curious behavior of the perturbation solutions is 
seen in the Lagrangian coordinates, not in the Eulerian coordinates. In order to have a 
more precise discussion, we have to transform the solutions into the form in the Eulerian 
coordinates. We will do so in a one- dimensional model in Sec. 0. 



IV. SECOND-ORDER SOLUTIONS 



In our previous paper, we derived the second-order solutions only for the case 7 = 4/3. 
In this section, we obtain the second-order solutions for the case 7 > 4/3 in an approximate 
form. To second order, Eqs. (0) and yield 



Vo x + 

a 



(!) (-A 1 ) 



+ 2-s 



(i) 
k,i) > 



(34) 



1 dP 

a 2 dp 



^1 + 2-4? - -^(Pb)v*</ - .s-v : i .si ■ : 2-.; 



2 (2) (l) / . a .(1) 



3 V l >3 



a 



1,3 



+ 1^(0^) ( s m v 2 s (1) + s (1) +s (1) V 2 s (1) +2s (1) s (1) ") 
a 2 dp v ' ,ij 9 j * ' jk j ' ik ' j q j '* ,j k Mj ' 

+ 4^(p b )p b w:]v 2 s « +s ^ (i) 



J^V-.l^-^S) 2 -^ 1 ^) , (35) 



where (•) $ denotes d/dqi. As in the first-order solutions, we decompose s^ 2 ^ into the longitu- 
dinal and the transverse modes as = V 9 C + C T with V q • £ T = 0. Then these equations 
are rewritten as 



V q x C T + 2-C T 

a 



-^-^(pb)eijkS,ijS7 2 q S M , 



(36) 



V 2 K + 2^-4 7 rGp b C-^(p b )V 2 ^ 



27iGp h 

1 d 2 P 

a 2 dp 2 



^■-j^-(Pb) (Vq^jV 2 ^, + SjjkS^jk + 2.V , ; V 2 .S' ,/ 



(p b )p b (v 2 SV 2 V 2 5 + V 2 ^V 2 5,, 



(37) 
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where we have neglected the first-order transverse perturbation S T for simplicity, and used 
Eq. (pUp . Taking the rotation of Eq. flSSP , we obtain 

- vi (c + 2^c) = ^(Pb) (s, ijk v 2 q s, jk + vjvj^ - vj^v;^ - s >jk v 2 q s, ijk ) . 

(38) 

The Fourier transform of Eqs. (|37D and (^) gives 

-|K| 2 fe+2^-47rG Pb C + — fa) \K\ 2 (] 
\ a a z dp J 



(2tt) 3 J-oo 

1 dP 



d 3 K'S(K',t) S(K-K',t) 



2nGp h {{K' ■ (K - K')f - \K'\ 2 \K - K'\ 2 ) 



+ ^^(p b ) {\Kf\K - Kf (K' ■ (K - K')) + (K' ■ (K - K')f 
+2\K - K'\ 2 (K f ■ (K- K')f) 

' <,rP (p b ) {\K'\ 2 \K - K'\ 4 + \Kf\K - Kf (K' (K — K' 



a 2 dp 2 



(39) 



a 



(2tt) 3 a 2 dp ^ W 7-oo 



d 3 K'S(K',t) S(K-K',t) 



K - K'\ 2 (K' ■ (K - K')) 



K[ IK' ■ (K - K') + \K- K'\ 2 } 



■{Ki - Kl) {K 1 ■ (K - K') + \K'\ 2 } 



(40) 



Using the Green functions G(K, t, t') and G T (t, t'), Eqs. Q3"9"D and (^D|) are solved in the form 

1 ft 



C(K,t) = 



\K\ 2 
1 



K 



dt' G(K,t,t') Q(K, t') 
dt' G T (t,t') 0f(K,t') 



(41) 



(42) 



where Q(K,t) and Qj(K,t) denote the right-hand side of Eqs. (|39|) and (|40"|) , respectively. 

In order to obtain an explicit form of the second-order solutions, we assume the Einstein- 
de Sitter background with a normalization so that a(t) = t 2 ^ 3 , and the equation of state as 
P = up 1 . The first-order solutions are then 



-1/3 



S oc const. , t 



S(K, t) oc r J- ±5/(8 _ 67) U\K\t-' +A l 3 ) for 7 ^ 4/3 



§(K,t)x t-V*±y/v>/M-B\K\* for7 = 4/3 



(43) 

(44) 
(45) 
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where 



.4 



1 



2C*2 g _ 2C*2 



|4 - 37 1 V Ci ' ~ 9d 
These first-order solutions yield the Green functions in the following form: 



G(K,t,t') 



71 



G T (t,t') = 3(t' - r 1 / 3 *' 4 / 3 ) , (46) 
J„ v (A\K\r^ 3 )J v (A\K\t'-^ 3 ) 



o , - 7 +-| t-^t'^ 
2 sin z/7T \ 3 



-J v (A\K\t~' y+4/3 )J- u (A\K\t'~' y+4/3 ) 



for 7^4/3. 



(47) 



G(K,t,t') 



1/25 „, -,,2^ X ^ --1/6+/7/6 /^- v /25/36-B|/<:| 2 ^/ v /25/36-B|/<:| 2 



2 V36 



. tv / 2 5/36-B|K| 2 t /-v/25/36-B|K|A for 7 = 



(48) 



where we have assumed that v = 5/(8 — 67) is not an integer. If we write the first-order 
solution as S(K,t) = D + (K,t)C + (K) + D~ (K ,t)C~ (K) , where D ± (K 1 t) are given by 
the form of Eqs. ([H]) and fl^5|), we obtain 
% 1 



C T (^,t) 



(2tt) 3 |Jf I 



d 3 K' E (K, K t)[C + (K')C + (K - K') 



hC + (K')C-(K - K') + C-(K')C + (K - K') + C~(K')C~{K - K' 
\K - K'\ 2 (K' ■ (K - K')) K[ IK' ■ (K - K') + K'\K - K'\ 2 ) 

-(Ki - K[) [K 1 ■ (K - K') + \K'\ 2 } 

d 3 K' (C + (K')C + (K — K') + C + (K')C~ (K — K') 



(49) 



(2tt) 3 \K\ 2 J-c 

+C-{K')C + {K - K') + C~(K')C~{K - K') 
■ E(K,K',t){(K' ■ (K - K')f -\K'\ 2 \K - K'\ 2 } 
+F 1 (K, K',t) {\K'\ 2 \K - K'\ 2 K K K') + (K' ■ (K - K')) 3 

+2\K - K'\ 2 (K' ■ (K-K')) 2 } 
+F 2 (K,K',t) \\Kf\K - K'\ 4 + \K'\ 2 \K - K'\ 2 K' (K-K' 



(50) 



where the time-dependent factors are given as 

f* dt' dP 



jW)) G T (t, (D+(K', t')D+(K - K', t' 



a 2 (t>) dp 

+D + (K', t')D-(K - K', t') + D~(K', t')D + (K - K', t' 
+D-(K',t')D-(K - K',t' 



(51) 
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E(K,K',t) = f ' dt'2-KG Ph {t')G{K,t,t')(D + (K',t')D + (K - K',t') 

+D + (K', t')D~(K - K', t') + D-(K', t')D + (K - K\ if) 
+D-(K', t')D~(K - K\ If)) , (52) 



F 1 (K,K',t) = f jL^ {ph{t '))G(K,t,t') (D + (K',t')D + (K - K',t>) 
+D + (K', t')D~(K - K'. if) + D~(K', t')D + (K - K\ if) 
+D~(K', t')D~{K - K', t')) , (53) 



F 2 (K,K',t) = f ^^{p h {if)) Ph {t')G{K^t'){D+{K\t')D+{K-K\t' 



-D + (K', i')D~{K - K', t') + D~(K', t')D + (K - K', t ; ) 
-D-(K',lf)D-(K - K',t')) = (7 - l^K, K',t) . (54) 



It is cumbersome to perform the integration of Eqs. (pTD-( p^) in a complete form unless 
7 = 4/3. (See Ref. [[18] for 7 = 4/3.) However, we can obtain the temporal factors in an 
approximate form in the following way. By the definition of the Bessel function, 

*** ( = ± ^TWTT) (T) • (») 

and thus if A\K\t~" /+4: ^ 3 <C 1, we can utilize the following approximation formulae: 

j ±v (A\K\r^) ~ { 4^TTT t±5/6 ■ (56) 

Note that these formulae are useful in the case 7 > 4/3, because they give the leading term 
with respect to t when 7 > 4/3. Substituting these formulae into Eqs. (pTj)-(|54|), we have 

(57) 

E?(K,K',t) ~ 3r( _ 1/4)2 l-^i ') Ont-3) for 7 = 2, (58) 

. . 15tt 1 ( A 2 \K'\\K - K'\Y 4/ „ , . 

E(K,K',t)~ rr—, r — 1 U t 4/3 , 59 

v ' ' ; 28sinz/7r(4-37)r(z/ + l) 3 r(-z/+ 1) V 4 / 

Fl(JC k' i) ~ A2(4 ~ 37) / ^jKiiK-xi y 

U ' ' } - 6smi/7r(5-2 7 )(10-67)r(i/ + l) 3 r(-^ + l) \ 4 J 

for 7 ^ 5/2, 5/3 , (60) 
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12sin(57r/7)r(2/7) 3 r(12/7) V 4 j \ 5, 

for 7 = 5/2, (61) 

F x {K,K',t) ~ -— M 2 I K/ H K ~ /3 _ ln A for7 = 5/3 . (62) 



80vr \ 4 / V5 

Let us re-examine the relation between the perturbative solutions and the polytropic 
index, which seems curious in the first-order level, as we mentioned at the end of the previous 
section. In the second-order level, the ratio of E(K , K' ,t) to the other temporal factors 
(e.g. Fi(K, K', £)) can be taken as a measure of the pressure effect, because E(K, K', t) is 
of gravitational origin and the others are of pressure origin, and thus the ratio is similar to 
the Jeans wavenumber Kj oc a^ 7-4 -' 7 ' 2 in the first order. The ratio reads 

E [ K ^ K '^\ „ A- 2 t 2 ^ 3 ~ a 37 " 4 . (63) 
F 1 (K,K',t) C 2 1 ; 

This means that the curious tendency of the first-order solutions is, unfortunately, unchanged 
at second order, contrary to our expectation. This result may be a consequence of the 
perturbation scheme we adopt. See Sec. [V| for detailed discussion on this point. 

V. ILLUSTRATION IN SOME MODELS 

In this section, we illustrate the perturbation theory formulated in the previous sections 



with examples in one- and two-dimensional systems. In our previous paper ||18|| ) we computed 
the power spectra of density perturbations in a one-dimensional model for the case 7 = 4/3. 
Here we calculate the power spectra for the case 7 = 5/3, and discuss the differece of the 
power spectra for the variation of the polytropic index 7. It is of significance to compute 
and compare the power spectra in the Eulerian coordinates, because the evolution of density 
perturbations has to be discussed in the physical Eulerian coordinates, and it is non-trivial 
how a physical variable is rewritten due to the transformation between the Lagrangian 
and the Eulerian coordinates. Moreover, we present the density field in a two-dimensional 
model, and clarify how the pressure effect appears in a spatial pattern of the density field by 
comparison with the dust case. In this section, we assume the Einstein-de Sitter background 
with the scale factor a(t) = t 2 ^ 3 for simplicity. The power spectrum of density perturbations 
is defined as V(k,t) = (\5(k,t)\ 2 ) , where k is a wave vector associated with the Eulerian 
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coordinates x, 5 = (p — Pb)/Pb is the density contrast, and (•) denotes an ensemble average 
over the entire distribution. 



A. Power spectra in a one-dimensional model 

We calculate the power spectra of density perturbations in a one-dimensional model for 
the case 7 = 5/3. We did this in our previous paper |L8| for the case 7 = 4/3. Here we 



choose another value of 7 and see the difference of the results for the variation of 7. The 
first-order solution is then written as 

S(K, t) = D + {K, t)C + (K) + D~(K, t)C~{K) , (64) 

where if is a component of the direction of inhomogeneities in the Lagrangian wave vector 
K and 

D±(K, t) = r 1 / 6 J- T5/2 [A\K\t-^) . (65) 

The Jeans wavenumber is found to be Kj = ^/Qt 1 /^ /A from Eq. (|3^) . 

We consider how to determine C ± (K) from the initial conditions for an illustration. 
Here we set the initial density contrast Si n and the initial peculiar velocity v m so that they 
coincide with those given by the Zel'dovich approximation, which is the Lagrangian first- 
order approximation for a dust fluid. The Zel'dovich approximation in a one-dimensional 
system is written as 

x\ = qi + t 2/3 ^,i(gi) , x 2 = q2, x 3 = q 3 , (66) 

= i+^Vnte) -1 ' (67) 

where ^(qi) is an arbitrary spatial function, describing initial inhomogeneities. Then we 
have 

1 



i+*,ii(<zi; 

V in = (V ini 0, 0) = 



u (gi), (68) 
V 3 ^,^), 0, 0) = (^,i(<Zi), 0, 0) , (69) 



t=t in 



3 



.3 

where we define an initial time t{ n = 1. As for the case P = Kp 5//3 , the first-order solution 
gives 

= K 2 [j^ /2 (A\K\)C + (K) + J 5/2 (A\K\)C-(K))] , (70) 
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vT n (K) = IK 



+ |-ij 5/2 (A|i^|) + 1 j 5/2 (A|if|r^)| 1 

Comparing Eqs. (|68[) and (|70|), and Eqs. Q69|) and ([71]), we find 



(71) 



C + W = " (cob(A|X |) - sin(A|^|)] 5? n (K) , (72) 



C-(if) = -\/^ ( sin (^l) + _Loos(A|Jir|)j 8 in (K) . (73) 

The initial density perturbation 8- m (K) = \5i n (K)\ exp(z0j<-) is chosen so that (^(i^)! 2 oc 
|i^| n with the spectral index n — 0, ±1, and the phases 0^ are randomly distributed on the 
interval [0, 2tt]. We set the constant A so that the Jeans wavenumber Kj is 80 at the initial 
time, t = t; n , where a — 1. 

To compute the power spectra within the Lagrangian approximations, we have to take 
caution about the difference between the Lagrangian and the Eulerian wave vectors, K 
and k. The Lagrangian solutions are obtained in terms of K, while the power spectra are 
presented by using k. Thus we have to transform the Lagrangian solutions into the form 
in the Eulerian space. The way of the transformation is described in, e.g. subsection 4.3 of 



Ref. [18 



In Fig. |l|, we show the power spectra V(k,t) at a = 1000, where Iris a component of 
the direction of inhomogeneities in the Eulerian wave vector k, using the Eulerian linear 
theory and the Lagrangian first-order approximation. Instead of the power spectrum it- 
self, we present the 'transfer function', V(k, t)/V(k, t- m ), for convenience because it does not 
depend on the initial conditions in the Eulerian linear theory but does in the Lagrangian 
approximations generally. The spectra by the Lagrangian second-order approximation are 
not presented, because they are almost coincident with those by the first-order approxima- 
tion, as in the 7 = 4/3 case. Indeed the difference between the Lagrangian first-order and 
second-order approximations in the 7 = 4/3 case is less than 10 % at \k\ ^ 150, and that in 
the 7 = 5/3 case becomes still smaller, less than 1 % at \k\ ^ 150 within our illustrations. 
(See Sec. (V| for the reason.) 

In our previous paper, we compared the Eulerian linear theory and the Lagrangian ap- 
proximations in the 7 = 4/3 case, where the Jeans wavenumber kj is a constant. In this 
case, the Eulerian linear density perturbations with wavenumbers smaller than a constant 

15 



wavenumber always grow, while those with wavenumbers larger than that always decay 
with acoustic oscillation because of the constancy of the Jeans wavenumber. On the other 
hand, in the Lagrangian approximations, small-scale perturbations are developed by the 
nonlinear effect, and as a result, the difference between the Eulerian and the Lagrangian 
approximations becomes large especially at high-frequency region. (See Fig. 2 of Ref. [0.) 

Now we observe the results of the 7 = 5/3 case, Fig. [I]. In this case, the Jeans wavenumber 
depends on time, and it becomes about 2500 at a = 1000 whereas it is set as 80 at the initial 
time. This means that the Eulerian linear density perturbations with wavenumbers between 
80 and 2500 are initially oscillating, but become growing modes later. Actually we can see 
this tendency at high-frequency region in Fig. |]. As for the Lagrangian approximation, 
small-scale perturbations are enhanced because of the nonlinear effect, as in the 7 = 4/3 
case. The difference between the Eulerian and the Lagrangian approximations is, however, 
not so large because of the behavior of the Eulerian density perturbations mentioned above. 

For comparison of the 7 = 4/3 and 5/3 cases in the Lagrangian first-order approximation, 
we show in Fig. |] the transfer function for both the cases, using the same initial conditions. 
This figure tells us that the growth of density perturbations computed by the Lagrangian 
approximation is suppressed by the pressure more weakly in the 7 = 5/3 case. This implies 
that the curious behavior of the Lagrangian perturbation solutions is preserved even if we 
observe it in the Eulerian coordinates. 

B. Density field in a two-dimensional model 

Next we consider an illustration in a two-dimensional model. In a dust model, Buchert 
and Ehlers |]19| showed the density field with the Zel'dovich and the 'post-Zel'dovich' approx- 
imations. Following their illustrations, we present a realization of the density field (mapped 
by 128 2 particles) with our approximations, in order to see how the pressure effect appears 
in a spatial pattern. We set the initial conditions for the scalar function S(q,t) as 

S(q,t in ) =^EE ^2 I y2 [^{K iqi + K 2 q 2 + <P{K U K 2 ))} , (74) 

Kf + Kl^Q, lifi i2 = 0,1,-.., 5, 

where the phases (fi(K 1 , K 2 ) are random numbers between and 2n, and the amplitude fi is 
chosen so that /i = 3.0 x 10~ 3 . The periodic boundary condition is imposed. We consider 
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the cases in which the equation of state is given as P = ftp 4 / 3 and P = /tp 5 / 3 , assuming the 
Jeans wavenumber Kj ~ 8. 

Setting the initial conditions at a = 1, the time evolution of the density field is shown 
in Figs. |3| and |j. In the 7 = 4/3 case, the evolution obviously proceeds slowly because of 
the pressure effect. In Fig. |3|, shell crossings just arise in the dust case, while the evolution 
remains still quasi-nonlinear regime in the 7 = 4/3 case (\5\ < 1.0). In Fig. |], shell crossings 
are being formed in the 7 = 4/3 case, while in the dust case high-density structures are being 
dissolved. In these figures, the difference between the first- and second-order approximations 
seems still small on large scales (compare (a) and (b), and (c) and (d)), although the second- 
order solutions should compensate shortcomings of the first-order approximation on small 
scales, as was discussed by Buchert and Ehlers Jl9| for the dust case. 

Above we have mentioned the 7 = 4/3 case, but what will happen if we take larger value 
of 7 such as 5/3? To answer this question, we show in Figs. || (e) and £| (e) the results 
computed by the Lagrangian first-order approximation in the 7 = 5/3 case. (The results 
by the Lagrangian second-order approximation are omitted, because we may presume them 



easily from other results presented.) As we stated in Sees. Ill and IV, the pressure effect in 
this case becomes weaker than in the 7 = 4/3 case. Indeed we see that the spatial density 
pattern resembles that in the dust case, rather than that in the 7 = 4/3 case. 

In our perturbation scheme, shell crossings arise in the 7 = 4/3 and 5/3 cases in spite of 
the presence of the pressure effect, but the features of first collapsing objects are manifestly 
different from those in the dust case. (Compare, e.g. Figs. |3] (a) and [| (c).) The growth of 
small-scale structures is particularly suppressed because of the pressure effect, and therefore 
the size of overdense region becomes larger, as if the density field was spatially coarse- 
grained. (Compare, e.g. Figs. ^ (a) and (c).) Consequently our perturbation scheme may 
work like the 'truncated Zel'dovich approximation' |20| [21], |22| , which yields a coarse-grained 
density field of the original Zel'dovich approximation. 



VI. DISCUSSION AND CONCLUDING REMARKS 

We have developed a perturbation theory in the Lagrangian hydrodynamics for a cos- 
mological fluid with pressure. Hydrodynamic equations in the Lagrangian coordinates have 
been solved perturbatively up to second order, extending our earlier work. In our earlier 
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work [fLaH , we solved the first-order perturbation equations in the Einstein-de Sitter back- 
ground, and the second-order ones explicitly for the case 7 = 4/3. In this paper, we have 
obtained the first-order solutions in non-flat backgrounds and flat backgrounds with A ^ 0, 
and the approximate second-order solutions for the case 7 > 4/3. We have found that in 
several cases, the first-order solutions are written in terms of Gauss' hypergeometric func- 
tion. We have also presented illustrations in one- and two-dimensional systems, showing 
how our approximation theory describes the evolution of cosmological inhomogeneities. 

In Sec. [V|, we have computed the power spectra of density perturbations in a one- 
dimensional model for the case 7 = 5/3 with the Eulerian linear theory and the Lagrangian 
first-order approximation, and have shown some amount of the difference between them. Our 
numerical calculations have also shown the difference between the Lagrangian first-order and 
second-order approximations, smaller than that in the 7 = 4/3 case. Let us investigate the 
reason of the smallness by considering single- wavemode perturbations and evaluating the 



ratio of the second- to the first-order solution, as we did in subsection 4.4 of Ref. [18|. We 



assume that the first-order solution is written as 

S (9ut) = -^Re [(c + (K)D + (K,t) + c~(K)D-(K,t))exp(iK qi )] , (75) 

where e is the amplitude of the initial density perturbations, c (K) denote constants of 
0(1), and D ± (K,t) are given by Eq. (H). Then, from Eq. (p0[), the second-order solution 
becomes 

e 2 

Cfa,*) ~ -—Re[F 1 {2K,K,t)exp(i2Kq 1 )) . (76) 

47T 

For a concrete estimation, we assume A\K\t-i+V 3 < 1 and use the approximation formulae, 
Eq. (|56|) , for the cases 7 > 4/3. This assumption is reasonable because this is equivalent to 
taking into account perturbation modes whose Lagrangian wavenumbers are smaller than 
the Jeans wavenumber. The first-order and second-order solutions are then reduced to 



V 



S( qi ,t) ~ 7^(4^J t 2/3 Re[exp(^ gi )], (77) 
/ A\K\\ 2u 

C(gi, t) ~ -e 2 A 2 4^ t" 2 ^ 4 Re[exp(z2K qi )] , (78) 



2 

where v = 5/(8 — 67), and thus we find 

C(?i,*) 



S( qi ,t) 



A\K\\ v+ \_^ nl , /K\ 2 fA\K\ 



< e (qfM] t - 2 7+ 10/3 „ £ ( ^_ ) t 2/3 (7g) 
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Note that the factor e(A\K\/2) v t 2 ' 3 corresponds to the Eulerian linear density perturbation 
and is of order unity at most in our case. Then we can show that \C/S\ <C 1, since the 
assumption A\K\t~ 7+i / 3 <C 1 is equivalent to \K\/Kj <C 1. 

In the above estimation, the second-order solution ((qi,t) is of purely pressure origin 
because of the one-dimensionality, and thus can be regarded as a measure of the 'second- 
order pressure effect'. Manifestly the effect of ({qi,t) becomes weaker in time as we take 
the larger value of 7. This curious fact is exactly the same as what we have addressed at 
the end of Sees. [TTT] and [TV]. Now let us examine the cause of the fact. We remark the terms 
of pressure origin in the perturbation equations, Eqs. (|20|), fl37|), and (pBp. Then we see 
that all the terms of pressure origin have time-dependent coefficients such as dP/dp(pt,) and 
d 2 P/dp 2 (pb)Pb ; which behave as 

dP d 2 P 

-fyiPb) K 7|pT(^)Pb oc a" 37+3 , 

under the assumption P oc p 7 . These coefficients originate from the perturbation scheme, 
and we can safely claim that these coefficients yield the curious behavior of the perturbation 
solutions. In addition, these coefficients will appear at any order in the perturbation scheme, 
and therefore the curious behavior will arise, i.e. the larger value of 7 will produce the weaker 
effect of pressure at any order, as far as we consider the Lagrangian perturbation scheme. Our 
two-dimensional illustration also indicates how the evolution of inhomogeneities is sensitive 
to the variation of 7; the pressure works effectively in the 7 = 4/3 case, but does not in the 
7 = 5/3 case, although it depends on the choice of values of parameters in general. Buchert 
et al. |13|] argued that the 7 = 2 case corresponds to the adhesion approximation ||10|| , but, 



considering our illustration, it seems difficult to realize the adhesion-like approximation in 
the 7 = 2 case within the Lagrangian perturbation scheme. 

However, there should be no such curious matter in the exact level of hydrodynamic 
equations. To see this, let us consider the one-dimensional case, where the relation between 
the Eulerian and the Lagrangian coordinates are given as 

Xi = ?i + Si(gi,t) , x 2 = q 2 , x 3 = q 3 . (80) 



Under the assumption P = np 1 ', the exact equation for si is [17, 18 



s\ + 2** - 4nG PhSl - K1P f ai -" = , (81) 
a a 2 (l + s^i)^ 7 
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where the fourth term of the left-hand side holds the pressure effect. This term also has 
the time-dependent coefficient, dP/dp(pb), but simultaneously includes the effect of inho- 
mogeneities by (1 + Si,i) 1+7 in the denominator. As long as | | <C 1, the results of the 
perturbation theory are reproduced, but once the flow lines of the fluid approach to the 
shell-crossing singularities, 1 + s^i — > 0, the effect of inhomogeneities becomes strong. In 
this situation, the larger value of 7 gives the stronger effect of pressure, and thus no curious 
matter will arise. 

In our earlier work, and also in this work, we have experienced the shell-crossing problem 
in spite of taking into account the pressure effect. However, we can expect that this problem 
will also be avoided in the exact level, because the fourth term of the left-hand side of 
Eq. (|8lD will become very large near shell crossing, 1 + s± 1 — > 0, and will stop the growth 



of density enhancement. (Some implication may be obtained by Gotz [23], who solved the 
one-dimensional exact equation for the case 7 = 1 without cosmic expansion.) 

The above discussion implies that we have to admit that our perturbation scheme yields 
some artificial results. This is true, but the Lagrangian perturbation scheme is a natural way 
to solve the hydrodynamic equations in cosmology, and our formulation will give a useful 
tool for large-scale structure formation in a practical sense. It is, in principle, applicable to 
any cosmological situation in which velocity dispersion arises and is written as a function 
of the density only. Actually Fig. |] has shown that our scheme works better than the 
Zel'dovich approximation beyond shell crossing, giving some kind of spatial coarse graining 



of the density field, as is given by the truncated Zel'dovich approximation |20], [2T|, |22 



Detailed analyses of comparison of our scheme and the truncated Zel'dovich approximation 
(and also the adhesion approximation) will be provided in a separate publication. 

As for the shell-crossing problem, Matarrese and Mohayaee |24| have treated it in La- 
grangian perturbative approach for two-component fluid. They also experienced shell cross- 
ing in usual perturbative Lagrangian approach, and introduced the 'stochastic adhesion' 
model to overcome the problem. It will be interesting to probe how to treat the dynamics 
when shell crossing is occurring, or how to avoid shell crossing by taking account of the 
pressure effect in a sophisticated manner. 
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FIG. 1: The 'transfer function' of density perturbations at a = 1000 computed by the Eulerian lin- 
ear theory, and Lagrangian first-order approximations. It does not depend on the initial conditions 
in the Eulerian linear theory, but does in the Lagrangian approximation. 
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FIG. 2: The 'transfer function' of density perturbations at a = 1000 computed by the Lagrangian 
first-order approximation in the 7 = 4/3 and 5/3 cases. Small-scale perturbations in the 7 = 5/3 
case are developed more effectively than in the 7 = 4/3 case. 



24 



dust 1st 



dust 2nd 



V 










'V 


V 






. 


■■■■ 


• 






■ 





(a) 

Y=4/3 1st 







.- 

\ 


\ 


■ 




* 

• 

i 




(b) 



y=4/3 2nd 





(c) 




Y=5/3 1st 


•■■ 






c 

\ 
















i 







(e) 

FIG. 3: The particular density field of a two-dimensional model at a = 1000. Shell crossings just 
occur in the dust case, (a) First-order approximation without pressure (the Zel'dovich approxima- 
tion), (b) Second-order approximation without pressure (the 'post-Zel'dovich' approximation), (c) 
First-order approximation with pressure, 7 = 4/3. (d) Second-order approximation with pressure, 
7 = 4/3. (e) First-order approximation with pressure, 7 = 5/3. 
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FIG. 4: The particular density field of a two-dimensional model at a = 3000. Shell crossings 
just occur in the 7 = 4/3 case, (a) First-order approximation without pressure (the Zel'dovich 
approximation), (b) Second-order approximation without pressure (the 'post-Zel'dovich' approx- 
imation), (c) First-order approximation with pressure, 7 = 4/3. (d) Second-order approximation 
with pressure, 7 = 4/3. (e) First-order approximation with pressure, 7 = 5/3. 
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